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The four-dimensional ±J random-bond Ising model is 
studied using ground-state calculations. System sizes up 
to N — 6 4 spins are considered. Here it is found that 
the ferromagnetic-spin glass transition occurs at a critical 
concentration p c = 0.28(1) of the antiferromagnetic bonds, 
which is comparable to values found previously by high- 
temperature series expansions. The transition is character- 
ized by a correlation-length exponent v — 1.0(1) and an order- 
parameter exponent = 0.4(1). Thus, this transition is in a 
different universality class from four-dimensional bond perco- 
lation, where v = 0.678(50) and (5 = 0.639(20). 

Keywords (PACS-codes): Spin glasses and other 
random models (75.10.Nr), Numerical simulation studies 
(75.40.Mg), General mathematical systems (02.10.Jf). 



INTRODUCTION 

The study of how order arises in nature is in the center 
of interest in several areas of physics, especially in ther- 
modynamics and statistical physics. A special focus is 
on the question, how order can emerge, even in or due 
to the presence of quenched disorder. One type of sys- 
tem, which has attracted a Jpt of attention during the 
last decades, are spin glassesO. 

In this work, the transition from one ordered phase, the 
ferromagnetic phase, to another ordered phase, the spin- 
glass phase, is studied for a four-dimensional random- 
bond Ising model. The model is studied by means of 
ground-state calculations, using the genetic cluster-exact 
approximation (CEA) method. This approach has the 
advantage that one does not encounter ergodicity prob- 
lems or critical slowing down like when using algorithms 
which are based on Monte-Carlo methods. To the au- 
thor's knowledge, there are no numerical studies of the 
four-dimensional random-Land Ising model (except at 
p = 0.5). However, for twcoQ and threeQ dimensions, the 
random-bond Ising model has been studied numerically 
at T — already, resulting in critical values p s c q — 0.10 for 
square lattices respectively p c c ub — 0.22 for cubic lattices. 

The model treated here consists oi N — L A Ising spins 
Cj = ±1 on a hypercubic lattice. The Hamiltonian of 
this system is given by 



H 



^2 JijViVj 

(i,3) 



(1) 



driven by increasing the fraction p of antiferromagnetic 
(Jij = — 1) (AF) bonds. For low concentration p, the 
system is ferromagnetically ordered, while for interme- 
diate values of p, spin-glass order arises. For very high 
concentration, the system shows AF ordering. 
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FIG. 1. Schematic phase diagram of the four-dimensional 
random-bond model. The diagram is symmetric with respect 
to p = 0.5, with an antiferromagnetically ordered phase for 
large p. SG denotes the spin glass phase. The multicritical 
point is denoted by (p*,T*). 

In previous work, the model has_been studied using 
high-temperature series expansionsUTJ. In Fig. |lj the 
schematic phase diagram is presented. The line connect- 
ing the multicritical point (p*,T*) with (p — p c ,T = 0) 
starts verticallyEHl9 at the multicritical point. Therefore, 
p c ~ p* can be expected. The values previously ob- 
tained for the concentration at the multicritical point are 
p* « 0.3 (Ref. |), p* = 0.28(1) (Ref. |) and p* = 0.290(1) 
(Ref. g). 

The rest of the paper is organized as follows: Next, the 
algorithm to generate spin-glass ground states in four di- 
mensions is briefly explained. Then, the results describ- 
ing the ferromagnetic-spin glass transition are presented 
and finally a summary is given. 



where (. . .) denotes a sum over pair of nearest neighbors. 
Systems with quenched disorder of the bonds Jij = ±1 
are studied. The ferromagnetic-spin glass transition is 



ALGORITHM 

From the computational point of view, the calculation 
of spin-glass ground states is very demanding, .-because it 
belongs to the class of the NP-hard problemsOii3. This 
means that only algorithms are available, for which the 
running time on a computer increases exponentially with 
the system size. Only for the special case of a planar 
system without magnetic field, e.g. a square lattice with 
periodic boundary conditions in at most one direction, 
exist efficient polynomial-time "matching" algorithmsEj. 

For higher dimensions., one must rely on exact methods 
like branch-and-bouncoll3 or branch-and-cutE3ii3, which 
are able to treat only small systems. For that reason, 
recently some heuristics-Jiave been introduced. By using 
a hierarchical approacblia one can calculate true ground 
states in four dimensions up to size L = 5, while for L p. 6 
it is not clear whether trua-sround states were foundtj. 

The method applied herec9, is able to calculate ground- 



states in four dimensions up to size L — 7 






tech- 



eaetic algorithmEJE3 and on 
3 which is an optimization 



nique is based on a special £ 

cluster-exact approximatior 

method designed especially for spin glasses. Now a brief 

description of the techniques are given. 

Genetic algorithms are biologically motivated. An op- 
timal solution is found by treating many instances of the 
problem in parallel, keeping only better instances and 
replacing bad ones by new ones (survival of the fittest) . 
The genetic algorithm starts with an initial population of 
Mi randomly initialized spin configurations (= individu- 
als) , which are linearly arranged in a ring. Then v x Mi 
times two neighbors from the population are taken (called 
parents) and tw(\-gffspring are created using the so called 
triadic crossoverO. Then a mutation with a rate of p m 
is applied to each offspring, i.e. a fraction p rn of the 
spins is reversed. Next, for both offspring the energy 
is reduced by applying the CEA. This algorithm con- 
structs iteratively and randomly a non-frustrated clus- 
ter of spins, whereas spins with many unsatisfied bonds 
are more likely to be added to the cluster. The non- 
cluster spins act like local magnetic fields on the cluster 
spins. For the spins of the cluster, an energetic mini- 
mum state can be calculated in. polynomial time by using 
graph-theoretical methodsLj~E3: an equivalent network is 
constructedEa, the maximum flow is calculatecOEII and 
the spins of the cluster are set to the orientations lead- 
ing to a minimum in energy. This minimization step is 
performed n mm times for each offspring. 

Afterwards each offspring is compared with one of its 
parents. The pairs are chosen in the way that the sum 
of the phenotypic differences between them is minimal. 
The phenotypic difference is defined here as the number 
of spins where the two configurations differ. Each parent 
is replaced if its energy is not lower (i.e. better) than the 
corresponding offspring. 

After this creation of offspring has been performed v x 
Mi times, the population is halved: From each pair of 



neighbors the configuration which has the higher energy 
is eliminated. If not more than 4 individuals remain, the 
process is stopped and the best individual is taken as 
result of the calculation. 

The whole algorithm is performed hr times and all 
final configurations which exhibit the lowest energy are 
stored, resulting in n g statistical independent ground- 
state configurations, i— 3y comparison with an exact 
method, it was shownEil that genetic CEA, with an ap- 
propriate choice of the parameters (Mi, v, n m ; n , p m ), 
indeed calculates true ground states. 

The probability that a certain ground-state configura- 
tion isJound by this method is not equal for all ground 
statesE3, i.e. the algorithm imposes a bias. In this work, 
the magnetization is the main quantity of interest. To 
test, how large the influence of this bias is, for L = 6, 
p = 0.27, where fluctuations of the magnetization are the 
largest, 64 realizations were considered. For each realiza- 
tion, (biased) ground states were generated using 200 in- 
dependent runs. Next, for each realization a (thermody- 
namically correct) set of ground states was generated in 
such a way that each configuration contributes with the 
same probability!^. For both sets of ground states the 
magnetization was evaluated. For the biased set an aver- 
age value of to = 0.501(11) was obtained, while the cor- 
rect thermodynamic result is m — 0.504(12). This shows 
that the influence of the bias on the result of the mag- 
netization is small. Please note that that is in contrast 
to other quantities like the distribution P(q) of overlaps: 
The overlap q a @ between two independent ground states 
{erf}, {of} is given by q = ^ J^i a t u i- When measuring 
the fraction aj(0.5) = J_ Q 5 P(q)dq of small overlaps, one 
obtains an average value x(0.5) = 0.051(7) for the biased 
set of states, while the from the correct thermodynamic 
treatment one obtains quite a different result 0.014(8). 

Hence, when restricting the measurement to quantities 
depending solely on the magnetization, it is sufficient to 
use the biased data. This allows much shorter running 
times, because only few ground states per realization are 
needed. 



RESULTS 

With-the genetic cluster-exact approximation, it it 
possibleEj to obtain ground-states for four-dimensional 
systems up to L = 7. Since the study of the transi- 
tion involves the calculation of a Binder cumulant, it is 
necessary to average over many samples of the disorder. 
For that reason, only sizes L < 6 were considered in this 
work. A mutation rate of p m = 0.05 was used. The other 
simulation parameters are shown in Tab. I. 

Systems with concentrations of the antifcrromagnetic 
bonds in the range p € [0, 0.32] were treated. For each 
realization tir — 5 independent runs were performed. 
The resulting ground-state energy as a function of the 
parameter p is shown in Fig. for different system sizes. 
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TABLE I. Simulation parameters: L = system size, Mi = 
initial size of population, v — average number of offspring per 
configuration, n m in = number of CEA minimization steps per 
offspring, r = typical computer time per ground state on a 
80MHz PPC601, N L = total number of realizations of the 
random variables. 



For small concentrations p, the ground state is mainly 
ferromagnetic. It follows that all the AF bonds are not 
satisfied, so the ground state energy increases linearly 
with p like e(p) ps — 4 + 8p. For larger concentrations, the 
ground-state energy approaches the p — 0.5 limit. With 
increasing L the ground state energy decreases, because 
the periodic boundary conditions impose less constraints 
on the system. For L — * oo and p — 0.5 ^ground state 
energy of eoo(0.5) = —2.095(1) was foundEJ. 




for the magnetization M = -^ ^ <Xj is evaluated. The 
average (...) denotes the average over different ground 
states of a realization, while [. . .] j is the average over the 
disorder. In Fig. y the Binder cumulant is shown for L = 
3,4,5,6. To keep the figure clear, only the largest error 
bar is shown. All curves intersect near p c = 0.28 ± 0.01. 
Because of the huge numerical effort, the simulations are 
restricted to small system sizes, which prevents a more 
accurate result. 

Only few studies concerning the four-dimensional 
random-bond model have been performed before. Us- 
ing high-temperature series expansions for the location of 
the multicritical point (p*,T*), values of p* ps 0.3 (Ref. 
|), p* = 0.28(1) (Ref. |) and p* = 0.290(1) (Ref. |) 
were determined. These results compare well with the 
value of p c obtained here, indicating that the the line 
connecting X&-, 0) with (j>*,T*) is vertical, or nearly so, 
as expectecoEEl. 
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FIG. 3. Binder cumulant g(p, L) of the ground-state mag- 
netization as a function of the AF-bond concentration p for 
system sizes L — 3, 4, 5, 6. For clarity, only the largest error 
bar occurring in the data is shown. Lines are guides for the 
eyes only. 



FIG. 2. Ground-state energy e(p,L) as a function of the 
AF-bond concentration p for system sizes L = 3,4, 5, 6. Error 
bars are much smaller than system sizes. Lines are guides for 
the eyes only. 

From Fig. |2| is clear that the energy as function of con- 
centration is not a suitable quantity for determining the 
critical concentration p c , where the ferromagnetic be 
ior disappears. For this purpose the Binder cumulant] 



g(p, L ) = 2 



3- 



KM 2 )} 2 , 



(2) 



For the Binder cumulant the following finite-size scal- 
ing relation is assumedo 



g(p,L)=g(L 1 ^(p-p c )) 



(3) 



By plotting g(p, L) against L x l v {j) ~ p c ) with correct pa- 
rameters Pc, v the datapoints for different system sizes 
should collapse onto a single curve near (p— p c ) = 0. The 
best results were obtained for p c — 0.28 and 1/v — 1.0. 
In Fig. [| the resulting scaling plot is shown. It is pos- 
sible to change the value of v in a wide range without 
large effects on the scaling plot. Thus, here a value 



of v = 1.0(2) is estimated. This is compatible with 
the value v =_0.8(1) found by a high-temperature se- 
ries expansion^ at the multicritical point. Although 
there is no direct correspondence between the T = 
ferromagnetic-spin glass transition and the finite tem- 
perature (p = 0.5) paramagnetic-spin glass transition 
at T c — 2.03(3), it is remarkable that by Monte-Carlo 
simulationsEj of four-dimensional spin glasses near T c 
a similar value of v = 1.0(1) was found. Finally, it 
should be pointed out, that the result given above is 
not compatible with the corresponding exponent for four- 
dimensional bond percolationtll of v = 0.678(50). 
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4. Scaling plot of the Binder cumulant as a function 
Vc)L 1/v with p c = 0.28 and v = 1.0. The line is a 



FIG 
of (p- 
guide for the eyes only. 



In Fig. the behavior of the average magnetization 
to = [(A/)] j is shown as a function of p for different 
system sizes_This quantity has the standard finite-size 
scaling forrr£3 



ro(p, L) = L- p l»fn{L l l u {p - p c )) 



(4) 



By plotting L^^ v ni{p, L) against L x l v {p—p c ) with correct 
parameters p c ,[3,v the datapoints for different system 
sizes should collapse onto a single curve near (p— p c ) = 0. 
The best result was obtained using p c = 0.28, 1/v =1.0 
and (3/v = 0.4. It is shown in Fig. g for L = 3,4,5,6. 
From variations of the value /3/v, the value of the expo- 
nent j3 = 0.4(1) is estimated. Again, this values signify 
icantly differs from the exponent (3 = 0.693(20) foundEZl 
for the order parameter of the four-dimensional bond- 
percolation. 




FIG. 5. Average ground-state magnetization m(p, L) as 
a function of the AF-bond concentration p for system sizes 
L = 3,4,5,6. All error bars are much smaller than symbol 
size. Lines are guides for the eyes only. 



SUMMARY 

The ferromagnetic-spin glass transition of the four- 
dimensional random-bond Ising model was studied at 
T = using ground-state calculations. The genetic 
cluster-exact approximation method was applied. Be- 
cause of the high computational effort, only small sys- 
tems of size N < 6 4 could be studied. Since the ground- 
state problem is NP-hard, it is very unlikely that signifi- 
cantly larger sizes can be studied in the near future. 

By evaluating the Binder cumulant, a critical concen- 
tration of p c = 0.28(1) was found, which is comparable 
to values which were obtained by high-temperature series 
expansions. 

Using a finite-size scaling analysis, the critical expo- 
nents for the divergence of the correlation length v = 
1.0(1) and for the magnetization (3 — 0.4(1) were deter- 
mined. Hence, the transition is clearly in a different uni- 
versality class than four-dimensional bond percolation. 
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